Magnetically Hidden Order of Kramers Doublets in d 1 Systems: Sr 2 VC>4 
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We formulate and study an effective Hamiltonian for low-energy Kramers doublets of d -ions on a 
square lattice. We find that the system exhibits a magnetically hidden order in which the expectation 
values of the local spin and orbital moments both vanish. The order parameter responsible for a 
time-reversal symmetry breaking has a composite nature and is a spin-orbital analog of a magnetic 
octupole. We argue that such a hidden order is realized in the layered perovskite Sr2VC>4. 
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The ions with odd number of electrons have ground 
states which are at least doubly degenerate. The lat- 
ter, known as Kramers degeneracy, is protected by time- 
reversal symmetry [lj|. In Mott insulators the exchange 
interactions between Kramers doublets often induce a 
time-reversal symmetry breaking phase lifting the local 
degeneracy. As a result, a magnetic state with a well 
defined ordered pattern of local magnetic moments is 
formed at low temperatures, which, in turn, leads to the 
new Bragg peaks in neutron scattering experiments. 

The above conventional picture may, however, fail for 
the insulating systems in which the local and non-local 
interactions between magnetically active degrees of free- 
dom compete with each other. This competition may 
lead to a spontaneous time-reversal symmetry breaking 
through the development of a more complex order pa- 
rameter, while the magnetic dipole moments of spin and 
orbital origin, being locally entangled, would remain dis- 
ordered across the transition. The resulting order can be 
thus hidden to some experimental probes. 

In this Letter we show that such a magnetically hidden 
order may be realized in transition-metal compounds. 
In the following we focus on a system of d 1 ions with 
partially filled t2 g levels on a square lattice, like the 
layered insulating compound Sr 2 V04 @, H, Here, 
a square lattice of V 4+ ions is formed in the a6-plane 
by corner-shared VC-6 octahedra, elongated along the c- 
axis. Sr 2 V04 undergoes a phase transition on cooling 
at around T c ~ 100 K. The crystal structure remains 
tetragonal across the transition, the c/a ratio jumps to a 
somewhat higher value, and the magnetic susceptibility 
shows a sharp drop near T c 0] • It has been suggested that 
a (weakly) first-order transition to an antiferromagnetic 
and orbitally ordered state could be responsible for the 
observed anomalies However, the neutron scattering 
experiments were not able to detect any magnetic order 
in the low temperature phase Q. The nature of the or- 
dered phase still remains experimentally unknown. The- 
oretically, as a possible candidate, a nontrivial orbital- 
stripe order coexisting with collinear antiferromagnetic 
spin order with large unit cell has been proposed [5(. 

Here, based on the microscopic theory, we propose that 
an unusual symmetry-breaking phase, induced by spin- 



orbit coupling and compatible with tetragonal crystal 
symmetry, is realized in Sr 2 V04. We argue that the stag- 
gered ordering of composite spin-orbital objects, mag- 
netic octupoles, is responsible for the phase transition. 
There is no static order of either spin or orbital mag- 
netic moments in the ground state, hence the absence of 
magnetic Bragg peaks. 

The octupolar order has thus far been discussed pri- 
marily in the context of /-electron systems 0, 0] ■ The pe- 
culiar physics of ti g orbitals in tetragonal compounds un- 
covered here shows that this unusual state is well hosted 
by d-electron systems, too. 

Low-energy effective Hamiltonian.- We first introduce 
the local degrees of freedom and then discuss the ex- 
change interactions between them. The V 4+ ion has a 
single unpaired electron residing in the ti g manifold of 
xy, xz, and yz orbitals. The tetragonal elongation of 
the oxygen octahedra along the z \\ c-axis only partly 
lifts the threefold orbital degeneracy: The xy orbital is 
pushed to a higher energy, while xz and yz orbitals re- 
main degenerate. The orbital angular momentum is un- 
quenched and the spin-orbit coupling is active. We thus 
start with a local Hamiltonian H$ = A c f(| — l 2 z ) — \l ■ S 
consisting of a tetragonal crystal field A c f and a spin- 
orbit coupling A. Here S is an electron spin, and I = 1 
is an effective angular momentum with \l z — Q) =\xy), 
\l z =±X) = — A=(i\xz) ± \yz)). The total magnetic mo- 
ment M = 2S + L, where L = —td is a true angular 
momentum and k is a so-called covalency factor of or- 
der one [l| . The eigenstates of Hq are spanned by three 
sets of Kramers doublets. The corresponding local level 
structure is schematically shown in Fig. []Ja). Concern- 
ing the energy scales involved, the ab initio study of 
Sr 2 VO"4 electronic structure suggests A cf ~ 80 meV [5], 
and A ~ 30 meV for free V 4+ ion is known experimentally 
1]. In what follows, we retain only the two low-energy 
levels and neglect the one located at high energy ~ A c f. 
We label Kramers partners within a doublet by isospin 
index f and J, (s z = ±4), while the two doublets are de- 
noted by pseudoorbital index ± (o~ z — ±1). The wave 
functions of the ground state doublet, a z = +1, are 
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FIG. 1: (Color online) (a) A schematic view of the local level 
structure: the t2 g level is split into three sets of Kramers 
doublets by a tetragonal crystal field A c f and a spin-orbit 
coupling A. (b) Density profile of an electron in the ground 
state doublet. The chirality of the isospin up wave function, 
associated with the gradient of its phase, is shown by small 
arrows. The electron spin S and orbital angular momentum 
L are coupled antiparallel to each other. 



The corresponding electron density profile has an axial 
shape, and the two Kramers partners, apart from the 
direction of an electron spin, are distinguished by a chi- 
rality associated with the gradient of their phases deter- 
mined by orbital angular momentum [see Fig.[]Jb)]. 

The first excited level, located at an energy <5 — A + 
i(2A c £ — A)(l/ cos 26 + 1), has the following eigenstates: 



|f)_ = sin0|-l,T) + cos0|O,i) , 
= sm0| + 14) + cos0|O,T) , 



(2) 



where tan26> = 2V2A/(A - 2A cf ). In the limit A cf > A 
of interest here, 26 ~ ir and 6 ~ A. 

We now discuss the interactions between neighboring 
ions. The effective Hamiltonian is obtained by projecting 
the corresponding t2 g superexchange model of Ref. [8[ 
onto the reduced Hilbert space spanned by the two low 
energy Kramers doublets ((T|) and We first consider 
the dominant part of the Hamiltonian and discuss later 
the effects of a finite Hund's coupling. We find 

H = J"£& ■ + ± af)(l ± aj) - | £V? , (3) 

(ij) 

where the first term describes the exchange coupling be- 
tween the neighboring states, the SU(2) isospin degrees 
of freedom are represented by Si, and Pauli matrices 
<7j denote pseudo-orbitals, referred to as simply orbitals 
from now on (not to be confused with original t2 g or- 
bitals). The upper (lower) sign is taken for a bond along 
a(b) direction, J — t 2 /U, where t is a transfer integral 
and U stands for the Coulomb repulsion. The local level 
splitting 8 between two Kramers doublets is given by the 
second term. 

The ground state properties - The interaction between 
the isospins Si in Eq. ((3]) depends on the orbital occupa- 
tions through the operator ( 1 ± of ) ( 1 ± crj ) . Since this op- 
erator is non-negative, the isospin exchange is purely an- 
tiferromagnetic (or equal to zero). This suggests that in 



the classical limit, (siSj) = —\, the expectation value of 
the first term in the Hamiltonian vanishes and all orbital 
configurations are degenerate. This extensive classical 
degeneracy, inherent to the coupled spin-orbital systems 
8], is lifted here by the second term of the Hamiltonian 
9( selecting uniform orbital order with a z = +1. The an- 
tiferromagnetic ordering of isospins Si is then stabilized. 

While the ground state in terms of staggered order of 
isospins may, at a first glance, look conventional, its phys- 
ical properties are in fact very unusual and depend on 
the spatial orientation of the order parameter m = (2sq) 
[Q = (tt,it) is the ordering wave vector]. To illustrate 
this, we express the isospin s in terms of the original, 
physical spin S and angular momentum I, operators: 



s z 



(4) 



The in-plane components of isospins are represented by 
a composite object, which is, remarkably enough, a spin- 
orbital analog of a magnetic octupole. On the other 
hand, the axial component is equivalent to the physi- 
cal spin. Thus, in-plane ordering of isospins (m _L z) 
corresponds to a staggered order of magnetic octupoles, 
while the axial one (m || z) corresponds to a magnetic 
dipolar order. In the former case, the local expectation 
values of both spin and angular moments are exactly zero: 
(Si) — (k) =0. This can be explicitly verified using the 
wave functions |T]) of the ground state doublet: the in- 
plane g-factor of this doublet g a b = 0. In the case of 
a dipolar order, m || z, we find again a somewhat un- 
usual picture: the value of ordered magnetic moments is 
strongly suppressed, (M z ) = 1- k<1, because of com- 
pensation between spin and orbital magnetic moments. 
In the absence of covalency, i.e. k = 1, both (M) and 
g c = 2(1 — k) vanish. 

The effective Hamiltonian ([3]), with isospin rotational 
symmetry, cannot select a direction of the order param- 
eter m: The latter can be rotated from a purely dipo- 
lar character to a purely octupolar one at no energy 
cost. In other words, a Goldstone mode describing the 
out-of-plane rotation of isospins corresponds physically 
to fluctuations between dipolar and octupolar orderings. 
However, the SU (2) isospin symmetry is only approxi- 
mate and is broken by Hund's coupling Jh, neglected so 
far. Finite Jh induces the anisotropy term, which in the 
present case is of easy-plane form 7i a (i,j) = — aJs z s z , 
with a = 2Jh/U « 1. It confines the isospins in the 
plane and selects the staggered octupolar order, with 
vanishing dipolar moments (Si) and (Li) on every site. 
The emergence of this highly unusual state in an appar- 
ently simple d 1 Mott insulator like Sr2V04 is surpris- 
ing, given that its single- hole counterpart, d 9 perovskite 
La2Cu04, is a conventional antiferromagnet. The origin 
of new physics here is due to an unquenched spin-orbit 
coupling, the significance of which is bein g recogn ized in 
the context of various phenomena 2, 11,12, 13|, 14 1 . 
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Excitation spectra.- We now turn to the excitation 
spectrum above the octupolar ordered state. We first 
study the Hamiltonian © at a mean-field level, decou- 
ple it into isospin and orbital sectors, and discuss later 
the interactions between them. We employ isospin (or- 
bital) wave representation for s, («?$) operators in terms 
of Holstein-Primakoff bosons 6, (ai) and diagonalize the 
harmonic part of the mean-field Hamiltonian. 

The isospin (intradoublet) and orbital (inter- 
doublet) excitation energies are given by Wk = 
2J s ^(l- 7k + a 7k )(l+7 k ) and ft k = y/8{8 + 8 J o7k ), 
respectively. Here 7 k = ^(cosfc^ + cosk y ), 
J s = J (I + crfa*), and J a = J(s % s 3 + ±). At 
T = 0, we estimate J s ~ J and J ~ — 0.08J. The 
dispersion relations are plotted in Fig. [2] (left) for the 
realistic values of the parameters (S = 2.2 and a = 0.2 
in units of J). The in-plane isospin excitations are 
gapless, while the out-of-plane gap (~ y/a) at T point is 
induced by easy-plane anisotropy protecting octupolar 
order. The weakly dispersive interdoublet excitations 
are centered around 6. 

The above elementary excitations correspond, in fact, 
to the fluctuations of rather unconventional degrees of 
freedom. The in-plane isospin waves represent octupo- 
lar excitations, not directly coupled to the conventional 
spectroscopic probes such as neutrons. To make some 
predictions for inelastic neutron scattering experiments, 
we now discuss the magnetic excitations. To this end, we 
express the local magnetic moment M$ in terms of isospin 
and orbital wave operators (bi and a%, respectively): 

Mf ~ (6, + b\){oi + a\) , Mf ~ (a 4 + a\) , 

M? ~ i(l-«)e* <lR «(& i -&t) . (5 ) 

Here, only the leading order terms in boson operators 
are kept, and the isospin order parameter along y-axis is 
assumed. All three components of Mi consist of fluctu- 
ating parts only. Hence, no magnetic Bragg peaks will 
be seen in elastic scattering. The inelastic response of 
the ordered state is rather nontrivial. The fluctuations 
of M v component, which have the largest spectral weight 
of order one, are given by orbital excitations only. We 
thus predict well-defined interdoublet excitations at en- 
ergies ~ 5 [upper curve in Fig. [2] (left)] to be observed by 
polarized neutrons. The out-of-plane magnetic response, 
M z component in Eq. j5]), corresponds to gapped out- 
of-plane isospin excitations. However, the small spectral 
weight ~ (1 — k) 2 would probably make it hard to de- 
tect them. The in-plane M x response is given by the 
composite isospin-orbital excitations [see Eq. flSJ)]. The 
conventional spin wave excitations, inherent to magneti- 
cally ordered systems, are replaced here by a two-particle 
continuum [the shaded area in Fig. [2] (right)]. The latter 
is bounded from below by O k +Q, the interdoublet exci- 
tation energy at momentum shifted by Q. 

The interaction between isospin and orbital waves, ne- 
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FIG. 2: (Color online) Left: The dispersions of elementary 
excitations along the direction X(ir, 0) — > M(n,n) — > V in 
the Brillouin zone. The lower (upper) curve corresponds to 
the isospin (orbital) waves. Right: The spectrum of in-plane, 
M x component, magnetic excitations. It consists of a contin- 
uum of composite, isospin-orbital excitations (shaded area), 
and a quasiparticle part (lower curve, the width scales with 
intensity). The corresponding diagrams are also shown: The 
solid (wavy) lines denote the isospin (orbital) excitations and 
the open dots stand for the three-particle vertex in Eq. (B). 
Energies are given in units of J. 

glected so far, dynamically induce a linear coupling be- 
tween in-plane magnetic and isospin excitations. The 
three particle processes, responsible for this dynamic hy- 
bridization, have the following form (in units of J): 

W = Ei ? 'p 6 q fek + ^( b -A + &k&- q )}(a p + al p ) , (6) 

where the matrix element rj p = (cosp x — cosp y ), and 
p = q— k. The second order perturbative correction from 
TC to the in-plane magnetic susceptibility \xx is given by 
the lower diagram in Fig. [2] (right). When this correction 
is included, the spectral function of in-plane magnetic 
excitations \xx receives also a quasiparticle contribution. 
It lies below the continuum and is given by 

Its intensity is peaked at A-point and dies out when ap- 
proaching to the ordering wave vector [see Fig. [5] (right)]. 
The Goldstone mode is thus invisible, reflecting the ab- 
sence of magnetic Bragg peaks. Along (1,1) direction, 
the hybridization matrix element vanishes and the quasi- 
particle peak is absent. 

Finite temperature transition - We now turn to finite 
temperatures and show that interplay between interdou- 
blet and intradoublet excitations leads to a first order 
phase transition. The in-plane ordering of isospins breaks 
their rotational symmetry around z-axis. Note that the 
uniform order of low-energy doublets, a = (of) = +1, 
conserves the tetragonal crystal symmetry, as xz and yz 
orbitals are equally occupied on every site. The scale of 
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Mott insulator. The present theory suggests a nontrivial 
magnetic excitation spectrum in the ordered state, that 
can be verified by neutron scattering experiments. The 
spin-resolved circularly polarized photoemission experi- 
ment would be another test of the present scenario. This 
technique measures the I • S scalar product [17j, which 
we predict to be ~ 0.5 in the ground state. Finally, we 
suggest that another candidate to exhibit octupolar or- 
der is Sr2NbC>4, in which a more pronounced spin-orbit 
couplin g is expected. Sr2Nb04 is known to be a Mott in- 
sulator [18j ], however, its low temperature magnetic prop- 
erties have not yet been reported. 



FIG. 3: The isospin m and orbital a order parameters, isospin 
exchange J s , and orbital gap A vs temperature. The broken 
lines mark corresponding jumps at the first order transition. 
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the transition temperature is set by the isospin exchange 
energy J s , which depends on orbital correlations. On the 
other hand, the orbital correlations themselves are deter- 
mined by isospin ones controlling the orbital exchange 
J a . The nonlinear feedback effects between two degrees 
of freedom can drive a first order transition. 

To substantiate the above picture, we employ a mean- 
field approach. The isospin order parameter m is given 
by a self-consistent equation m = tanh(m J s /T) with 
J s = J (I + afcTj), where the orbital correlation func- 
tion (of of) = Ek7kV[^ktanh(fi k /2T)]. The energy of 
orbital excitations Slk is controlled by orbital exchange 
J = J(l - m 2 )/4. 

Shown in Fig. [3] are the temperature dependences 
of isospin and orbital order parameters, together with 
isospin exchange and orbital excitation gap A. The first 
order jumps are clearly seen in all physical quantities. 
This is compatible with the abrupt changes in lattice 
parameters and magnetic susceptibility observed experi- 
mentally near T c ~ 100 K [HEB]. 

Fig. [3] shows that not only isospin order parameter m 
but also the exchange energy J s jumps at the transition. 
This suggests pronounced magnetoelastic effects and ex- 
plains an abrupt in-plane contraction of the crystal, en- 
hancement of a c/a ratio. Note also that jumps in a and 
A imply an increase of the population of the low-energy 
doublet. As the latter does not include planar xy orbital, 
Jahn- Teller coupling acts in the same direction as mag- 
netoelastic effect and further enhances a c/a ratio. Fi- 
nally, since the ground state doublet |T]) is non-magnetic, 
the magnetic susceptibility drops at the transition to the 
level determined by the Van Vleck contribution from the 
transitions to the high energy doublets [l6| ■ 

In conclusion, a magnetically hidden octupolar order 
may be induced by spin-orbit coupling in d 1 transition 
metal oxides. We have argued that such a hidden order 
is realized in the perovskite Sr2V04, explaining thereby 
the puzzling absence of magnetic Bragg peaks in this 
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